Spectral Properties of Atoms in Fields: A Semiclassical Analysis 



P. N. Walker and T. S. Monteiro 
Department of Physics and Astronomy, University College London, Gower Street, London WCIE 6BT, U.K. 

(February 5, 2008) 



> 
o 



> 
oo 
O 
O 



o\ ■ 
"a : 

I ■ 

o : 



■ 

X 



We develop a semiclassical theory for the spectral rigidity of non-hydrogenic Rydberg atoms in 
electric fields and evaluate the significant deviations from the well-known Poissonian behaviour in 
the hydrogenic case. The resulting formula is shown to be in excellent agreement with the exact 
quantal result. We also investigate diamagnetic atoms; we find that, in contrast to the classically 
integrable atoms, diffraction has a small effect on the spectral rigidity in the classically chaotic atom. 
We show our predictions can also be of use in the mixed phase space regime. 
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The short-range spectral properties of non-hydrogenic 
Rydberg atoms in external fields were recently found to 
have an unexpected character [Q, the nearest neigh- 
bour distribution was neither Poisson nor Wigner-Dyson, 
but close to a new generic intermediate class known as 
half-Poisson. These findings have been investigated ex- 
perimentally for helium atoms in electric fields fol- 
lowing much interest in the dynamics of non-hydrogenic 
atoms in weak fields and the effects of core-induced chaos 
l^-Q. Concurrently, there has also been much inter- 
est in such intermediate Nearest Neighbour distributions 
with broad application in problems such as electron- 
electron interactions in closed mesoscopic devices [p|,p^, 
the metal-insulator transition pl] |, intruder states in nu- 
clear physics and inclusions in billiards ||]. These 
studies of intermediate statistics have exclusively em- 
ployed quantum calculations. 

Here we obtain spectral rigidities of non-hydrogenic 
Rydberg atoms in electric and magnetic fields from an 
accurate quantal calculation. We find interesting and 
substantial deviations from Poissonian (hydrogenic ) be- 
haviour in the integrable (electric field) and near inte- 
grable (weak magnetic field) case. We develop a semi- 
classical theory for Stark spectra which is in good agree- 
ment with the quantal results and shows that one-scatter 
diffractive orbits account for most of the effect. In con- 
trast we find a comparatively small effect in the case of 
fully chaotic Rydberg atoms. To our knowledge this is 
the first analysis of the curve- form of the spectral rigidity 
for a generic atom, which we derive from classical dynam- 
ical information. Our semiclassical analysis should be ex- 
tendable to the mixed phase space Kolmogorov-Arnold- 
Moser (KAM) case. 

The spectral rigidity, defined as 



A(L) = min ^ 
^ ' A,B L 
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where Af{E) is the spectral staircase function, p is the 
density of states (assumed constant over the range of L 
to be considered) and {■ . .)e denotes averaging over the 
spectrum, was first analysed semiclassically for classically 
integrable and chaotic systems by Berry JTst . It provides 
a measure of long-ranged deviations from the Weyl rule 
for spectral density, and for diffraction-free systems is 



controlled (for L > 0(1)) by classical dynamical correla- 
tions. For L ^ 0(1) it converges universally to the form 
A{L) — L/ 15 for symmetry reduced spectra. 

The hydrogen atom in a static electric field (of strength 
F) or in a magnetic field (of strength B) provided some 
of the cleanest illustrations of integrability, mixed phase 
space and chaos in a real system . The dynamics 
of the electron is two-dimensional and has a useful scal- 
ing property: the classical dynamics depends only on a 
scaled energy e — Ek^, where k = F^i/^ for the elec- 
tric field case and n = B~^/^ for the magnetic field case. 
This property is exploited in both experiment and the- 
ory. Spectra are obtained at fixed e and the correspond- 
ing eigenvalues Ki represent effective values of with 
fixed classical dynamics. Hydrogen in a magnetic field is 
near-integrable for e < —0.5; as the field is increased it 
makes a gradual transition to full chaos at e ~ —0.1. The 
Hamiltonian for hydrogen in an electric field is always 
separable. We consider field values where the eigenval- 
ues are well below the ionization threshold at e = — 2, so 
the system may be considered integrable and bound. 

Most experiments in fact investigate atoms other than 
hydrogen (typically He, Li or Rb). The useful scal- 
ing property may still be exploited but the inner core 
of electrons yields non-trivial effects: additional weak 
spectral modulations and spectral statistics near the 
Wigner-Dyson limit even for the integrable/near inte- 
grable regime fj^l^. Hence the interest in so-called 
core-induced chaos. The additional modulations are ac- 
curately described by Diffractive Periodic Orbit Theory 
[Q. An investigation ||l| of the NNS statistics for the 
lowest 40,000 eigenvalues showed that they are only near 
Wigner-Dyson for the lowest ~ 1000 states. For small h 
they were found to make a transition to an intermediate 
distribution P{s) ^ ase^"^^ with (a ~ 3 — 4), near to the 
Half-Poisson distribution (for which a = 4). 

Rydberg atoms and molecules are well described by 
Quantum Defect Theory, one of the most widely used 
theories in atomic physics . The core is described by 
a set of phaseshifts (quantum defects) 5i in each partial 
wave, quantifying the departure from pure Coulomb be- 
haviour. We consider the s-wave scattering case, where i5o 
is the only non-zero phaseshift, which describes lithium 
((5o ~ 0.4l7r) and helium {5q ~ O.Stt for triplet and 
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5q ~ 0.147r for singlet helium ) extremely well. 

The insertion of a single scattering channel can be rep- 
resented as a perturbation by a projection operator. In 
this case, it has been shown [QjIUJi^ that the eigenval- 
ues of the perturbed system remain trapped between the 
unperturbed eigenvalues. This trapping puts a strong re- 
striction on perturbations to long-range spectral correla- 
tions; indeed, in the limit h —^ , L oo, the perturbed 
spectral rigidity can differ from the unperturbed one by 
at most some value, bounded by and 2. 

We have calculated the lowest 36000 states of the Stark 
spectrum for magnetic quantum number m — 0, e = —3 
and a range of quantum defects. The results for (5o = 
(hydrogen) and Sq = 7r/2 ('lithium') for values of the 
effective k < 600 are shown in Fig ^ These val- 
ues span and go well beyond typical experimental val- 
ues. ( the experimental NNS statistics [|| correspond 
to K 130 - 150). Fig H shows that for our k < 600 
range, deviations from Poissonian/hydrogenic {6q — 0) 
behaviour are substantial. Even for such small values of 
the effective h though, it is evident that the perturbed 
curves are not obtained by a simple constant shift. 
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FIG. 1. Rigidities calculated from quantal spectra of hy- 
drogen and lithium-like atom for e = —3. The eflective values 
of hr^ (ie the k) are indicated. The straight line indicates the 
I//15 Poisson limit. The figure shows that the non-hydrogenic 
rigidities follow a significantly different curve. Surprisingly, 
given the h dependent nature of the diffractive process, the 
perturbed curves are essentially independent of fi below the 
saturation value of L ~ imax. The curves from our final 
semiclassical formula,Eq.(9), corresponding to both the SC 
(x = .16) and QM (x = .21) values of the constant shift are 
superposed on the k — 600 non-hydrogenic quantal rigidity, 
showing the excellent agreement and the relatively modest 
effect of the uncertainty in the shift. 

The electronic core yields a combination of Coulomb 
and s-wave scattering. The Coulomb scattering gener- 
ates Gutzwiller Periodic Orbits (POs): geometric orbits. 
The effect of the s-wave (Sq) scattering is to generate 
diffractive trajectories, 0{^/h) in amplitude weaker than 
isolated geometric POs. For atomic core-scattering or- 
bits which close at the nucleus corresponds to either peri- 
odic or half-periodic orbits : hence -for chaotic or regular 
dynamics- every diffractive orbit resulting from a single 
scattering is paired with a geometric periodic orbit or 
a half-periodic orbit of the same action. In the Stark 



case there are no half-period contributions; in the mag- 
netic field case there are half-period diffractive contribu- 
tions (the 'D' orbits seen in 0). All multiple scattering 
diffractive orbits can likewise be associated with a given 
combination of geometric periodic or half-periodic orbits. 
In contrast, in billiard systems diffractive orbits are gen- 
erally unrelated to the geometric periodic orbits. Their 
proliferation relative to geometric orbits in the chaotic 
regime is not restricted by this 'pairing'. It results in a 
non- vanishing semiclassical contribution [ p7[ . 

For integrable atomic spectra, the geometric contri- 
bution to the staircase function Mg{k) = ■^ci^) + 
Af^'^ln) includes a sum over contributions ||l8| from in- 
tegrable tori. In the scaled atomic spectra, AfQ'^{K) ~ 

\/«Ej(A/^j>'^^'''""'''^^^ The amplitudes Aj are the 
Berry- Tabor amplitudes for resonant tori which, in action 
angle variables is Aj = ^2Tr/\u}.dIj/dTj\Dct{de/dIj} 

All tori contain just one PO that collides with the nu- 
cleus. Trajectories that can diffract are therefore isolated 
in the usual sense. The amplitude Ad of the contribution 
of a particular single-scatter isolated diffractive trajec- 
tory of action S, in our 2D system ,was given in ||^ = 

S'j/7r-\/27r/mi2Ksin5oe*'^° sin0 where mi2 is an element 
of the reduced monodromy matrix M, (jj is the angle of in- 
cidence of the orbit relative to the field direction. In inte- 
grable systems toi2 (x S. Hence Ad scales as y/S/n. We 
investigate the ratio of the diffractive to corresponding 
geometric contribution, A-'jj/Aq = AdTj / {nAj) which 
we write as A^j^/Aq = iCj sin Soe"^^" Sj / k. By a unitary 
transform to action angle coordinates we can show |l|] 
that Cj ^ lu[^''ll!2 ^ where the uj are the frequencies of 
motion along the two independent degrees of freedom. 
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FIG. 2. Fourier transforms of scaled spectra for hydrogen 
and a lithium-like atom(5 — 0.5n) for e — —3 from a fully 
quantal calculation. The inset amplifies the region around 
S ~ 17. The amplitudes of the diffractive orbits may be 
obtained from the (complex) difference of the hydrogen and 
lithium traces. The amplitudes of multiple-scattering diffrac- 
tive orbits 0(h^^^) are relatively weak, but in the very large 
S regime (L ~ 0{1)) the proliferation in the number of n— 
scatter orbits with increasing action means that the n = 1 
contribution ceases to be predominant. 

Hence the Cj are independent of k and do not scale 
with orbit length. They are nearly constant and fluctu- 
ate weakly about an average value |l^ . From the Fourier 
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transforms of scaled spectra , we can confirm this be- 
haviour. The Ad/Ag ratios and hence the Cj statistics 
are obtained numerically from the first ^ 50 pairs. 

Since the ratio of the one-scatter diffractive orbits to 
the resonant tori contributions is 0{S/k), which in the 
unsealed spectra is 0{hT), on classical timescales diffrac- 
tion has a very small effect, but on quantum timescales 
(T ~ l/?i) both contributions are of the same order. We 
show that a semiclassical analysis of the rigidity is suffi- 
cient to reproduce such essentially quantum phenomena. 

We have also calculated the lowest 10000 states of the 
diamagnetic problem e — —0.6 (near-integrable) , and 
e = —0.1 (chaotic). The results are in Fig H] (inset). 
In the chaotic case we find almost no perturbation of the 
hydrogenic result. Isolated hyperbolically unstable orbits 
contribute like |2 — TrM|~^/^, which vanishes exponen- 
tially with the orbit period. In such a system, mi2 typi- 
cally diverges with the same Lyapunov exponent as Tr M. 
Hence A-'j^/A-'q is 0{^/fi) regardless of trajectory length. 
Proliferation of diffractive orbits relative to geometric 
ones is restricted to additional half-period contributions 
by some orbits, in contrast to the billiard result [ pT| , 
where any fraction may contribute. We conclude that, 
in the accessible h range, the diffractive effect remains 
small for the chaotic atom, but represents a substantial 
effect in the integrable case, on quantum timescales. 

In the KAM system, deviations qualitatively similar 
in form to the integrable case occur. We do not at- 
tempt a rigorous analysis of systems with a mixed clas- 
sical phase space, but the numerics support the notion 
that the structure of the Berry- Tabor formula ||l^ is re- 
tained for near-integrable KAM systems and that 
our analysis could also find application in such systems. 

The preceding arguments are of course incomplete for 
the rigorous development of a semiclassical limit, as we 
do not consider multiple scattering effects such as the 
creation of the combination orbits. However, we will see 
that in the numerically accessible regime (and indeed far 
beyond experimental resolution) it is not necessary to 
consider such effects to accurately reproduce the diffrac- 
tive corrections to the spectral statistics for L ^ C(l), 
up to a small constant correction. 

We now analyze the case of a non-hydrogenic Stark 
atom. Since we will only be interested in long orbits, 
for which the uniformity principle po| can be invoked to 
connect period with action, we can write the Berry- Tabor 
formula as a sum over all periodic orbits j (including 
retracing to negative time) in the form: 



AA°«=(«;) = sinJoT^-r-^ y (;.^.^^.g«(S,«-a,^/2±5o) (3) 
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where the phaseshift ±1^0 is positive for forward time trac- 
ings, and negative for negative time retracings. 

The rigidity formula (P, which is easily adapted for 
the problem in hand where the levels are the k^, leads 
to integrals involving the products (A/'gA/'g)k, {■N'gJ^d)k 
and {NdJ^d)k- These products involve many rapidly 
oscillating terms, which, following Berry [|l3|, we assume 
to vanish upon averaging. This approximation is usually 
known as the diagonal approximation, but here we retain 
different (off-diagonal) trajectories with identical action, 
which occur due to diffraction. 

The resulting equation for the rigidity can be written 



A*''(L,L,) = L/I5- 



where 
A^^(i,ic) 



(L,L,)-I-A^^(L,L,) , (4) 



-4.2sin^5o(#:) E;Q/o 
X dS{A]IS)5{S - Sj)G{l^^LS/2) (5) 



A^^(L,L,) = 2ns\n^5, C^^^'^' 

X dSA]S{S - Sj)G{A^LS/2) , (6) 

the sum is now only over positive traversals, and G{y) = 
1 - F^iy) - ^F'iy)? (Fiy) = siny/y), is Berry's orbit 
selection function, which is similar to the step function 
Q{y — tt). We have had to introduce an upper cut-off 
to the integrals, which will be seen to diverge. This 
divergence is a direct result of neglecting the rapidly 
oscillating terms, and will be discussed further below. 
We now concentrate on evaluating the integrands of (||) 
and (h). For long orbits, we invoke the Hannay-Ozorio 



de Almeida sum rule |20|, which can be expressed as 



\i-ais^ooJ2j ^"j^iS - Sj) = (A«/A2)/27rK3 which is in- 
dependent of K. We consider this to be accurate for 
5 > S**, so that our result for (Q) will only be accurate 
for L < Lmax = 2tt//S.hS*. Indeed Lmax also marks the 
onset of non-universal deviations of the non-diffractive 
result from the simple L/15 dependence Inserting 
the sum rule into (^ and (^, averaging over the distri- 
bution of Cj and evaluating the integrals, leads to closed 
form solutions that can be expressed in terms of special 
functions |l6| . Here, we only write out the asymptotic 
approximations to the solutions: 



A^«(i,L,) 



ttkAk 



[In. 



Lr. 



■ IE 



(7) 



Here Sj, aj, are the scaled action and the Maslov in- 
dex, Aj is the scaled amplitude factor from ||l^, A^^ = 
dAf'/dKU, and = dJ\f'/de\^. Note that A«/A, 

scales as k. 

The contribution to the spectral staircase from the cor- 
responding single-scatter diffractive periodic orbits yields 
a sum over the same set of periodic orbits as in Eq.(0): 



A^^(L,L,) 



2(C2)sin2(5o 



(8) 



where is Euler's constant. kA^ is a constant, indepen- 
dent of K and ^ —2e for our Stark spectra. The result 
is valid for L > Lc, and the ambiguity in Lc leads to 
an ambiguity in the constant term. For L <ti Lc both 
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and A-^*^ vanish. To proceed, one should invoke 
the semiclassical sum rule due to Berry . It is natural 
to then identify Lc as 2tt/A^Sc, {eg Lc ^ .94 for the case 
in Fig. ^) where Sc is the point where the semiclassi- 
cal and quantum asymptotes for the form factor coincide. 
Since we have neglected both off-diagonal corrections and 
higher order scattering contributions, we cannot expect 
to be able to evaluate the constant term correctly. In 
this case we can simply set Lc to unity, and accept that 
our result may not be accurate around L ~ 1. We note 
that in the GOE case it is not possible to evaluate the 
constant term semiclassically either p3| . 
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FIG. 3. Comparison between exact diagonalisation (QM) 
and Eq.(^) (SC) for the integrable case of atoms in electric 
fields for scaled field e = — 3 for k ~ 600, x = -21. In order to 
expose the approximately logarithmic nature of the diffractive 
correction here we plot the difference between the hydrogenic 
and difi'ractive rigidities. The inset shows the correspond- 
ing (quantal only) results for atoms in magnetic fields for 
the near-integrable/mixed case (solid line, e — —.6,k, ~ 150) 
which is qualitatively similar to the integrable case and the 
fully chaotic rigidity differences (dotted line e — — .1,k ~ 125) 
which are ~ 0. 



The final result for the rigidity is then 



A^°{L) = A°(L) 



sin^ 5q {Cj 



InL 



Ml 
2eL 



X (9) 



for Lniax ^ ^ 1, where x is the constant shift. For 
L<1, A^°{L) = L/15. The final formula hence does not 
contain the effective Planck's constant, proving the most 
surprising feature of Fig. |] namely that despite the h- 
dependent nature of the diffractive corrections, our quan- 
tal results (below £max) are essentially ?i-independent. 

We find (Cj) ~ 1.8, and (C^) ~ (Cj)^, for the e = 
—3 case. We stress that these were not free parameters. 
With Lc - 1. from @) and (|) we find x -16 for Sq = 
tt/2. However we estimate from the quantal results that 
X ^ 0.21. The correction required can be attributed 
to the neglect of both higher order scattering, and off- 
diagonal contributions. 

Our result is compared with the fully quantal values 
shown in Fig. ^ and the quantal and semiclassical rigidi- 
ties are plotted in Fig. . The agreement (to within 



the small correction to x ) is extremely good up to Lmax, 
where the breakdown was expected. The divergent na- 
ture of Eq. (||) for L < 1 is clearly seen. The sin^ Sq 
dependence has been verified by considering several dif- 
ferent defect values (not shown) |16|. 

To summarise, we have combined the semiclassical the- 
ory of diffraction and atomic Quantum Defect Theory 
with the Berry- Tabor trace formula to give diffractive 
corrections to the spectral rigidity of atoms in fields. We 
show that within a small constant shift, the semiclassical 
one-scatter results agree extremely accurately with the 
fully quantal results. 
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